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Abstract 



, Charge transport in tunnel junctions between singlet anisotropically paired 

| superconductors is investigated theoretically making use of the Eilenberger 

equations for the quasiclassical Green functions. For specularly reflecting 
tunnel barrier plane we have found and described analytically characteristic 
CN . singular points of the I-V curves which are specific for the case of anisotropic 

, pairing. All four terms for the electric current are examined. Two of them 

describe ac Josephson effect and two correspond to the quasiparticle current 
(the last term occurs only for a time-dependent voltage). Various momentum 



■ dependences of the order parameter in the bulk of superconductors and differ- 



ent crystal orientations are considered. The results of numerical calculations 
for the total I-V curves are presented for some particular cases. 



a 

I. INTRODUCTION 

Theoretical description of the Josephson and quasiparticle currents through a tunnel 
junction was developed many years ago for the case of s-wave isotropic superconductors. 
In particular, the total tunnel current under the externally applied time- dependent voltage 
was investigated in detail microscopically making use of tunneling Hamiltonian method |l],[2j 
(see also [Q and references therein). Now it becomes clear that in the case of anisotropi- 
cally paired superconductors the measurements of the Josephson effect and the quasiparticle 
current provide an important information about the structure of the superconducting order 
parameters from both sides of the junction and about the proximity and some other specific 
surface superconducting effects at the tunnel barrier plane j|-|J. The investigations in this 
field are of great interest and Josephson and quasiparticle tunneling experiments studying 
the anisotropic structure of superconducting order parameter for various high-temperature 
superconductors have attracted now much attention (see, for example, [|T0| ,|lT|j ) . 



One of the characteristic peculiarities of microscopic description of charge transport 
through tunnel junctions in anisotropically paired superconductors is, in fact, the failure 
of the tunneling Hamiltonian method applied to this case fl2| . This approach, being usu- 
ally quite suitable for tunnel junctions between s-wave isotropic superconductors, contains 
ambiguities in the case of anisotropic pairing due to substantial momentum dependence of 
the matrix elements describing tunneling between superconductors. One can show that the 
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choice of the tunneling matrix elements as being independent on the momentum direction 
(standard for the s-wave case) leads to confusing results for anisotropically paired supercon- 
ductors. Furthermore, an unambiguous choice of this dependence cannot be done within this 
method. This leads to the necessity of making use of a microscopic description of tunneling 
between such superconductors based on matching of the electron propagators at the tunnel 
barrier, since this approach leaves no space for ambiguity. 

Below we develop a microscopic study of the singular points on current-voltage charac- 
teristics for ac Josephson effect and for quasiparticle current of tunnel junctions between 
anisotropically paired superconductors under the externally applied ( and, generally speak- 
ing, time-dependent) voltage. In our microscopic approach we are based on the Eilenberger 
equations for quasiclassical electron propagators combined with corresponding boundary 
conditions and the microscopic expression for the tunnel current. In the particular case 
of s-wave isotropic superconductors our results coincide with those found by Larkin and 
Ovchinnikov [I]]. They obtained general expressions and described the singular points for 
functions I m (V) (m = 1,2,3,4), which enter the relation for the total tunnel current (see 
below, for example, Eqs.(^)-(0)). Here h^(V) describe the amplitudes of two terms in the 
expression for the Josephson current, which are reduced for time-independent voltage to 
Ji(V) sin (xi -X2 + 2eVt/h) + I 2 (V) cos (xi - X2 + 2eVt/h). Further, I 3A (V) describe the 
quasiparticle current. The function I^iV) occurs only in the case of time-dependent voltage 
while for permanent voltage the quasiparticle current is entirely reduced to the only term 
I 3 {V). It was shown, in particular, that at \eV\ = A\ + A2 the singular parts of functions 
h,4:{V) diverge logarithmically, while the functions I2,3{V) undergo a jump. Analogous singu- 
larities were found for time-dependent voltage oscillating with the frequency hu = Aj + A2. 



The singularity of I\ is known as Riedel singularity |13| . It is closely associated with the 
corresponding singularity in the density of states for superconductors at Twj = A. 

It is clear that Riedel's anomalies themselves must be washed out in the case of strongly 
anisotropic pairing when the density of states doesn't have the divergence. Nevertheless, as 
we show below, some new characteristic singularities appear even in the presence of nodes in 
the order parameter on the Fermi surface. The characteristic behaviors of the I-V curves turn 
out to be strongly dependent upon the crystal orientations of the superconductors, which, 
in particular, govern the spatial dependence of the order parameters near the barrier plane. 
But even in the simplest case of uniform distribution of the superconducting order parameter 
on both sides of a tunnel barrier the anomalies become strongly modified as compared to 
the isotropic case. The point is that for anisotropically paired singlet superconductors the 
functions I m (V) become also depending upon the momentum direction on the Fermi surface 
and enter the expression for the tunnel current as corresponding integrands. For obtaining 
the current-voltage relation in the vicinities of singular points one should carry out the 
integration over the momentum directions. We show that the nonanalytical behavior of 
the I-V curves takes place only for the voltages \eV\ which are equal to the values of the 
expressions ||A 2 (p 2 )| ± |Ai(pi)|| at their extremal points. Here pi is the direction of the 
incident momentum and p2 - of the transmitted momentum, which is directly connected 
with pi and the shapes of the Fermi surfaces. All the expressions are considered, for example, 
as functions of pi. Further, the singular behaviors of the I-V curves at these extrema turn 
out to be strongly dependent upon the type of the extremal points. This was pointed out 
earlier in where the characteristic behaviors of the I-V curves have been determined 
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for the important particular case of the quasiparticle current at low temperatures, for the 
crystalline orientations when there is no surface pair breaking. 

The situation becomes more complicated for the crystal orientations, for which the spa- 
tial dependence of the order parameter near the barrier plane plays an important role. The 
possibility for existence of quasiparticle bound states located in the vicinity of barrier plane 
is of importance in this case. At least two factors may be associated with such kind of effects. 
The former is the zero energy quasiparticle state which occurs at the barrier plane under 
certain general conditions PjMHiTjl . The latter is the appearance of additional quasiparticle 
bound surface states due to the spatial dependence of the order parameter, which is sup- 
pressed by the surface ||. We discuss below these possibilities and respective consequences 
for the current- volt age characteristics. 



II. MICROSCOPIC EXPRESSION FOR THE TUNNEL ELECTRIC CURRENT 

We consider a tunnel junction with transparency coefficient D < 1 and specularly 
reflecting barrier plane between two clean singlet superconductors. The external voltage 
V{t) = $2^) _ is assumed to be applied to the junction. Let the normal to the junc- 
tion barrier plane be directed along the x-axis: n || Ox. Then the microscopic expression for 
a tunnel current density in the first order of the barrier transparency D may be represented 
in the form f|T~ 



V X 1>0 J 



-tf-tf-tf'9?) (t,t) . (1) 

Here and below we put e = 1 as well as h,c = 1. 

The retarded, advanced and Keldysh quasiclassical matrix propagators are taken in 
Eq.(|TJ) at the junction barrier plane and must be calculated in zeroth order in junction 
transparency (i.e. for impenetrable barrier). They depend upon the respective momentum 
directions pi and p>2. Index 1 (2) labels the left (right) half space with respect to the bound- 
ary plane, and v x is the Fermi velocity component along the normal to the plane interface 
n. The integration in Eq.(P is carried out over the part of the Fermi surface with v x > 0. 
The relation between the incident and transmitted Fermi momenta ( that is between pi 
and P2 ) is as follows. The components parallel to the specular plane interface are equal to 
each other, while the values of normal components are determined by the condition, that 
Pi and P2 lie on respective Fermi surfaces. Naturally, in the particular case of identical 
superconductors with the spherical or cylindrical Fermi surfaces (provided the cylindrical 
axes are parallel to the barrier plane) the total incident and transmitted momenta are equal 
to each other p! = p 2 . 

Notations used in Eq.([l|) are in correspondence with the following example: 

t 

(g?-g?) (t,t)= J dt x ${t,t x )'${tx,t) . (2) 
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Here g RA ' K {t, = S t (t) g R ' A ' K {t - t x ) S+ (t x ) (Z = 1,2), where 

Si(t)=[ e e _^ )/2 J , Xi(t) = xi-2J ®i (f) dt' (3) 

and xi is the phase of the order parameter of Z-th superconductor at the boundary plane for 
the case of zero electric potential $j. 

Below we consider only singlet types of anisotropic pairing, for which matrix propagators 
may be represented in the form 

• (r ',) • w 

For further calculations it is also important that nonequilibrium effects, as a rule, are not 
essential in tunnel junctions. The voltage V simply shifts the Fermi levels in the electrodes 
relative each other by the value V. Besides, for finding the current in the first approxima- 
tion with regard to transparency, the calculation of the Green functions may be performed 
for impenetrable half spaces (disregarding the transmissions of electrons through the junc- 
tion). Under the conditions we consider the distribution functions of electrons remain the 
equilibrium ones and the effect of electric potential results only in the appearance of cor- 
responding spatial independent terms in phases of superconducting order parameters and 
Green functions on both banks of the junction (see, for instance, Eq.(|3|)). Then the following 
relationship is valid 

^H = (^»-^»)tanh(|,), (5) 

According to general symmetries of the propagators one has g A = f 3 (j) R ') T3, where 
f 3 is the third Pauli matrix. We suppose that the complex phase of the order pa- 
rameter within the superconducting half space with impenetrable boundary is constant 
in the coordinate and momentum spaces. Then this phase is just the quantity xi i n 
the Eq.([5[) and the Green functions further must be calculated for the real values of 
the order parameter, though its sign may depend on the momentum direction. Under 
these conditions one has f R ' A (— p, x, u) = —f +R ' A (p,x,uj), g R,A (—p,x,u) = g R,A (p,x,u), 
f R ' A (p, x, — uS) = f R ' A *(p,x,u), g R ' A (p, x, — u) = —g R ' A *(p,x,uj). We take into account as 
well that a propagator taken at the impenetrable boundary for the incoming momentum p 
is equal to the one taken for the outgoing momentum p, and v x (p) = —v x (p). Then one 
can get that for the time-dependent voltage V(t) = Vq + a cos(a;ot) the expression (fj) for 
the tunnel current acquires the form 



jx= Y\ J n ( — ) I ji (V + w„) sin ( xi ~ X2 + 2V t + — sin (u t) + u n t) + 



+j2(Vo + uj n ) cos ( xi - X2 + 2V t H sin (uj t) + to n t ) + 
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+j3(V r + u n ) cos (u n t - ^- sin (uJ tjj + j 4 (y + oj n ) sin (j n t - ^- sin (o; t)^| . (6) 
The following notations are introduced here 

UV)= J ^^(pOU^pO ,m = 1,2,3,4 , (7) 



Pi) = / ^ tanh Jm (/* ( W - 7) (tf* ( W ) + ft H) + (1 - 2)) , (8) 



oo 

W Pi) = tanh (^) ife (/* ( W - V) (/+* («) + /f («)) + (1 - 2)) , (9) 



CO 

h(V, Pi) = / ^ (tanh (^-) - tanh Qi)) Img? (u) ■ Img* (u - V) , (10) 



Wpi) = -/^T tanh (^)(ife^f (w-V) Img* (cu) + Reg^ (u-V) Img* (to)). (11) 

— oo 

The quantities Ii^iV) are even functions on V and 72,3 - odd functions. Note, that the 
interchanges / <-> / + in the expressions @, @ lead to the same result for the Josephson 
current after the integration over the momentum directions in (0). Eqs.(|6|)-(|TT|) are in 
correspondence with those obtained by Larkin and Ovchinnikov Q in consideration of s- 
wave isotropic superconductors. 

In the case of permanent voltage through the junction one should let a = in Eq.(^). 
Then one gets 

UV) = Ji(V) sin (xi -X2 + 2Vt) + j 2 (V) cos ( X i - X2 + 2Vt) + j 3 (V) . (12) 



It follows from Eqs. (|i2l) and (y), that the specific singular points on the voltage- 
current characteristic for the tunnel current for both time-dependent and permanent voltages 
through the junction are defined by the singular points of functions j m (V). If corresponding 
singular points are V — V m , then in the case of permanent voltage the specific points on 
the current-voltage characteristics are Vi 2,3 , while for the current © peculiarities appear 
for the values V m ,n = V m — uojq (m = 1, 2, 3, 4; n — 0, ±1, ±2, . . .). For large enough values 
n ((a/uo) < n) the amplitude of the current becomes quite small (according to respective 
behavior of the Bessel functions). 
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III. QUASICLASSICAL GREEN FUNCTIONS AT THE IMPENETRABLE 

BOUNDARY PLANE 



As it follows from previous section, for the calculation of the tunnel electric current one 
should consider the retarded electron propagators for the impenetrable superconducting half 
space and then find the values of the propagators at the boundary plane. The total analytical 
solution of this problem is quite complicated and hasn't yet been obtained for any particu- 
lar pairing potential leading to an anisotropically paired superconductivity ( excluding the 
particular orientation for which there is no surface pair breaking and the order parameter 
doesn't manifest a spatial dependence). The problem may be essentially simplificated if one 
is interested only in the singular points on the I-V curves and consequently the singular 
points of the propagators taken at the barrier plane. For the consideration of this problem 
we make use of the Eilenberger equations for the retarded quasiclassical propagators, which 
may be written for the case of superconductors with singlet pairing as follows (further we 
omit the superscript for the retarded propagators) 

' (2u + iv x d x )f(p, x, uj) + 2A(p, x)g(p, x, u) = 

< (2uj - iv x d x )f+(p,x,uj) - 2A*(p,x)g(p,x,u) = (13) 
. iv x d x g(p, x, u) - A(p, x)f + (p, x, uj) - A*(p, x)f(p, x,uj) = . 

For the sake of definiteness let us assume that the superconductor occupies a half space 
x > 0. 

Apart from a self-consistency equation for A(p, x), Eqs. (|l3"D have to be supplemented by 
a normalization condition 

g 2 + f-f + = -* 2 , (14) 

and boundary conditions for quasiclassical propagators. For the specularly reflecting impen- 
etrable surface one gets 

g(p, w) = g(p, u) \ x=0 , f(i>,u) = f(p,v)\x=o, f + (i>,u) = f + (p,uj) \ X=Q . (15) 

Here p is the direction of incident momentum and p - the direction of reflected momentum. 
The behaviors of the propagators in the depth of the superconductor 

g(p,u)\ x=00 = , / \ x=00 = -/ [p,U)\x=oa= I (16) 

must be also taken account of as the additional condition for the solutions of Eqs.(|T3|). 

Supposing that one can choose the gap function A to be real within the superconducting 
half space with impenetrable boundary (i.e. in the absence of a current across the junction), 
we define 

h = \(f-f + ) , f2 = Uf + f + ) ■ (17) 



2 v- * / > 2 
The Eqs. (|l3"l) , QUI) , being applied to /i >2 , take the form 
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{ 2ufi + iv x 8 x f 2 + 2Ag = 
d x g = -1^/2 , 



(18) 



g 2 + ft-fl 



-IT 



(19) 



The boundary conditions at x = for functions /i^ are the same as (|15D, while in the 
depth of the superconductor we have 



fi (p, 



Aoo(p)| 2 -^ 2 



2 a;=oo 







(20) 



The representations for g and / taken on the boundary, which turn out to be useful from 
the point of view of the consideration of singular parts of the propagators, may be derived 
from these equations. For obtaining the representations it is convenient to introduce the 
following function 



A(p,w) 



f™ A(p,x)f 2 (p,x,uj)dx 



jo°° f2(p,x,u)dx 

Since from the second and third equations of (|IBD one easily obtains 



(21) 



-0(0) = / A(p, x)f 2 (p,x,u)dx , /i(oo) - /i(0) = / f 2 (p,x,u)dx (22) 

V r J0 v r JO 



2iu r°° . 



the quantity A(p, oj) may be written also in the form 

0(oo) -g(0) 



A(p, to) = —uj 



/i(oo)-/i(0) 



(23) 



Here and below we denote 
g \ x =o= g(0) , g 



g{oo) , / | x . =0 = /(0) , / U=oo= /(oo) . 



After the substitution of the expressions Qlq) , (P0|) for the propagators in the depth of 
the superconductor into the Eq. ([2"3"D , it reduces to the following relation between ^(0), /i(0) 
and A(p, to): 



/m 7T A{p,uj)A 00 {p) - uj 2 A(p,u) 
0v u J - _ — t, = =^ : — /h u J 



UJ 



A^pW-uj 2 



uj 



(24) 



Entirely analogous relation may be written also for the momentum direction p, after 
that the boundary conditions (|T5|) allow us to write down the following representations for 
g(0) and /i(0) separately 
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7r A(p,w)A(p,w) 


A oc (p)A(p,u;) 


wA(p,u)-A(p,u) 


_A(pW|A 0O (p)| 2 -u; 2 



Aoo(p)A(p,^) 



or 



(25) 



/i(0) 



7T 



A(p,w) - A(p,u;) 
Eq. 



A 00 (p)A(p,u;) 



a; 



A QO (p)A(p,u;) 



a; 



A r 



P 



A oc (p)| 2 -^ 2 



(26) 



One can see from Eg. (|25|) that candidates for the singular points of the propagator 
g(0) are uj = 0, ±|Aoo(p)|, ±|Aoo(p)|, and, generally speaking, the singularities of functions 
A(p, uj), A(p, uj). Analogously, one finds from Eq.(^) that candidates for the singular points 
of /i(0) are uj = ±|Aoo(p)|, ilA^p)!. As far as the solutions of the equality A(p,u;) = 
A(p, uj) are concerned, we note that the consideration of Eq.fl2~3|) for the momentum direction 
p doesn't result then in the independent relation as compared to Eq. fl24|) for the momentum 
direction p. Thus, some additional information is needed for the consideration of these 
limiting cases in Eqs.fl2"5|), fl2"6|). 

Let us consider firstly singular parts of the propagators g s , f s , /+ taken at the boundary 
in the vicinity of the point uj = 0. Taking the low-frequency limit uj — > in Eqs.(|25|), (p6|), 
we find that the function /i(0) has no singularity at uj = 0, opposite to the propagator 
g(0) which turns out to have the pole at this point ( under the condition that A^p) and 
Aoo(p) have opposite signs): 



9.(0) 



1 7rA(p,0)A(p,0) 
wA(p,0)-A(p,0) 



(s^n(A 00 (p)) - s^n(A t>0 (p))) 



B B (P) 



UJ 



(27) 



It is remarkable that in addition to this relation one may find also the explicit expression 
for the quantity A(p, 0) through the inhomogeneous distribution of the order parameter. 
Indeed, since the function fi has no singularity at uj = (see Eq.(|l8|)), we can derive from 

±^/2s (when uj — * 0) between the singular 
) we get the equation for 



the normalization condition (19) the relation g s 
parts of g and ^2- Then from the last equation of the system 
the singular part of f 2 : d x f 2 s ± — A s = 0. The solution of this equation which satisfies the 
conditions (f20|), (|15|), is 



/ 2)S (p, x, uj) = / 2)S (p, 0, uj) exp 



2s5fn(Aoo(p)) 



A(p,x')dx' 



(28) 



Substituting this solution into Eq. (|2l|) we obtain 



A(p,0) = - 





v x 


\sgn( Aoo( 


P)) 




fexpf 2S ^( A : (P)) , 
Jo \ \v x \ J 


f* A(p,x')dx' 


dx 



(29) 



Note, that the determination of sign in Eq. (^) allows to fix sign in the relation between 
the zero frequency singular parts of the quasiclassical propagators in the close vicinity of 
the boundary: 



f s (p,x,u) = f+(p,x,w) = -i sgn(v x A^(p))g s (p,x } uj) . 



(30) 
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Thus, Eqs.fl2?D, ( p9|) and ( j30|) provide the quite general description for the zero frequency 
singular parts of the propagators taken at the boundary. If one is interested in only the 
singular points of the I-V curve, one can separate the problem of the solution of the self- 
consistency equation for particular pairing potentials for other investigations. The latter 
problem being an important part of the total theoretical description of the I-V curve for the 
tunnel junction, is very cumbersome and obviously includes large numerical investigations 
within the framework of any microscopic model for the pairing potential. 

Passing to the other candidates for the singular points uj = ±|A oc (p)|, ±|Aoo(p)|, we 
find firstly one important relationship for the function A(p, u;) at this frequency. For this 
purpose the asymptotic behavior of the function f 2 (p, uj) at x — » oo is of interest, which 
follows from Eqs.©-©: 



J2 (p, x, uj) oc exp I - j— j x \ . (31) 

One can see from Eqs.(|31|), ( pT|) that in the limit |a>| — > |Aoo(p)| the main contribution 
to the integrals in fl2ll) comes from the depth of the superconductor, where A(p, x) is equal 
to its bulk value. Then one obtains the relationship 

A(p,u/) -> A oc (p) for H->|Aoo(p)|. (32) 

Now one can consider the behaviors of g and fi in the vicinity of the points \u\ — > 
|Aoo(p)|, |Aoo(p)|. Taking the limit |a>| — > |Aoo(p)| or |Aoo(p)| in Eqs.(p5|), (||), one can 
find with the help of (133) the cancellation of divergences and, in the first approximation, 



the appearance of the square root nonanalitycal behavior of the form J A^ (p) — uj 2 (or 



A^, (p) — uj 2 ) for g(0) and /i(0) (and therefore /(0)) at these points. The exclusions 

may represent the orientations for which A(p) = A(p), when strictly speaking Eqs.(p5|), 
(PB|) themselves do not provide more information than (E^). It is well known that if this 
condition holds for all momentum orientations, there is no surface pair breaking and the 
divergences in the propagators at |a>| — ► |Aoo(p)| occur in this particular case: 

Z , ; * m ■ (33) 

v/|A(p)|' 2 -^ 2 v/|A(p)| 2 - w 2 

At last one should discuss the possibility for existence of quasiparticle bound states 
with nonzero energy located near the boundary. They may appear, for instance, due to 
the spatial dependence of the order parameter which is suppressed at the boundary. The 
bound state may be interpreted as a bound state in the "potential well" formed by the 
order parameter [§]. Since the quasiparticle bound state corresponds to the pole in the 
quasiclassical propagators, we simply add below the pole-like term to the singular parts of 
the propagators. If the propagators <?(0), /i(0) have similar poles describing localized states 
with nonzero energy, then it follows from Eq.fl2"3p and the boundary conditions (|T5D, that at 
the point of the pole A(p, uj) = A(p,u). 
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Taking into account the results obtained above, the nonanalytical terms of g(0) and /(0) 
may be written in the form 



Q 9 (p) 



uj + iS uj — h(p) ■ sgn(uj) + iS 



+ C (p, u ) sgn(u) JA^ (p) - cu 2 + C (p, u) sgn(uj) J A^ (p) - cu 2 + 



(34) 



Is (P,V) x=0— — J" H T7TT 7 s 

c<j + 20 uj — n,(pj • sgn(uj) + to 

+ E (p, W ) ^ (p) - W 2 + E (p, w) v/A^ (p) - W 2 + . . • , 5^+0 . (35) 

The following relationships take place: Bf(p) = — s^n(v a; A 00 (p))_B 9 (p) = —Bf(—p), 
Q g (p) = \Q f (p)\, Q f (p) = -Q}(-p), h(-p) = h(p). 

Below we are interested only in the values of functions Q g j{p) near the poles uj = ±h(p). 
In considering positive and negative poles one should take into account that quantities 
Q g (p,u), ReQf(p,u) are even functions, while ImQf(p,uj) is an odd function of uj. We do 
not indicate further the explicit energy dependence of functions Q g j, using only their values 
at the positive pole. 

Since for \uj\ < rain (|Aoo(p)|, |Aoo(p)|) the quasiparticle density of states in the contin- 
uum is zero for fixed momentum direction, it is natural to demand Ira C, Ira Q g = under 
this condition. In the particular case A(p, x) = —A(p,x) , when A(x = 0), /i(0) = 0, one 
may obtain from the Eilenberger equation Ref(0) = for the frequencies < |Aoo(p)|. 
Then for this frequency interval one has Re (E (p, uj) + E (p, uj)) — 0. 

As we show below, the location and type of singular points on the I-V characteristics 
of the tunnel current associated with the quasiparticle bound states, are determined, in 
particular, by the extremal and nonanalitical points of the function h(p). From the boundary 
conditions for the propagators and the parity of pairing we get h(p) = h(—p). Taking into 
account also that for p || n it follows p = — p, we find that the function h(p) (as well 
as the total propagator g(0)) must have an extremal value for the direction p along the 
normal to the boundary n. Analogously, we get that the function Qf(p) takes purely 
imaginary value for the direction p || n. In particular, it follows from here for the the 
orientations A(p, x) = —A(p,x), that Qf, g (p) = for the momentum direction p || n. 
The other characteristic points of the function h(p) are the momentum directions where the 
quasiparticle state bound to the boundary disappears. 

The results of numerical calculations of quantities B g ((f>), Q g {4>) = \Qf(<f>)\, ReQf(<p) 
and h{4>) are presented in Fig. 1. We consider a tetragonal superconductor with cylindrical 
Fermi surface and impenetrable specularly reflecting interface at x = (the cylindrical 
axis z is parallel to the boundary plane). The pairing interaction is supposed to have 
the particular form V ((/>,(/)') = 2Vcos(20 — 20 o ) cos(20' — 20 o ), which results in a d-wave 
order parameter. Here is the azimuthal angle in xy - plane, which is counted from the 
direction of the normal to the boundary. Angle <p describes orientation of the crystalline 
rro-axis with respect to the normal to the boundary. Then the order parameter has the form: 
A(0, x) = A(x) cos(2(/> — 20 o ), where A(x) has to be evaluated self-consistently. The angle 
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defines the direction of incoming momentum along the quasiparticle trajectory. We choose 
O = 7r/9, T = 0.45T C , A /(2T) = 2, where A = A(x = oo), T c is the critical temperature. 
For the surface to lattice orientation O = 7r/9 the order parameter is suppressed in the 
vicinity of the boundary up to the value A(x = 0) = 0.28A . Quantities B g ((f>), Q g ((/)), 
ReQf((f>) are normalized in Fig. 1 to the value vrAo, while the quantity h(<ft) - to the value 
Ao- One can see that the bound states with nonzero energy ±h(4>) exist only in the narrow 
angular region <p G (—0.095, 0.095) in the vicinity of the normal to the boundary. The mid- 
gap states exist within two other wider regions of momentum directions, where B g (<j)) ^ 0. 
It follows from Fig. 1, that the maximal value of the function h((f>) is h m = 0.7A and the 
value of h(4>) at the edge points is h e d = 0.63A . It is worth noting that for the momentum 
directions p (<f> e d = ±0.095), where the bound states with nonzero quasiparticle energy 
disappear, the following relation holds h e d = min (\A 00 (p)\,\A 00 (p)\). This means that 
discrete bound levels combine with the continuum spectrum at the momentum directions <p e d- 
Quantities C (p, uj), E (p, uj) (or C (p, uj), E (p, uj) - depending on what quantity is smaller 
among |Aoo(p)|, |Aoo(p)| ), which are present in Eqs.(j34p, (|35|), diverge at the momentum 
directions 4>ed and the energy \uj\ = h e d- 



IV. CURRENT- VOLTAGE CHARACTERISTICS FOR THE JOSEPHSON AND 

QUASIPARTICLE CURRENTS 



The peculiarities of the I-V curves are due to the singular points of functions I m (V, pi). 
If these functions don't depend upon the momentum directions, their singular behavior 
with the variable V directly describes the respective behavior of the I-V curves, according 
to Eqs. (H) - (PP . In contrast to this case, for anisotropically paired superconductors 
the integration over momentum directions greatly modifies the singularities of the current- 
voltage characteristics. Then the singular points of functions j m (V) describe the singularities 
on the I-V curves. The other distinctive feature of the anisotropically paired superconductors 
which essentially influences the behavior of the Josephson effect as well as the quasiparticle 
current, is their sensitivity to the inhomogeneities and interfaces. The suppression of the 
anisotropic order parameter at the boundary results in the specific surface superconducting 
effects. At least two important features must be taken into account in this context. The 
former is associated with the possibility for the opposite signs of the bulk order parameter 
taken for the directions of the incident and reflected momenta. This results in the zero energy 
quasiparticle bound states at the specular boundary. The latter appears if the additional 
quasiparticle bound states localized near the barrier plane occur due to the particular form 
of the spatial dependence of the order parameter (see, |T9|,p|,|20|,p^-p!7[] ) . 



A. Crystal Orientations with no Surface Pair Breaking 



Specific features of current-voltage characteristics for the tunnel electric current for 
anisotropically paired superconductors differ from the s-wave isotropic ones even in the 
case of conventional boundary conditions, when the surface doesn't suppress the supercon- 
ducting order parameter. This holds if the values of the order parameter taken for the 
incident and reflected momenta are equal to each other for all momentum directions. The 
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consideration of this particular case is just the subject of current section. Let us consider 
the behavior of functions j m (V), which determine the current- voltage curves, for the cases 
when the electron propagators at the impenetrable plane have the form fl3"3|) coinciding with 
their bulk expressions. Substituting these propagators into Eqs.(|^)-([TT|), we get 

-\u-V\) 
— ■ — x 



J 1 (^P 1 ) = -A 1 (P 1 )A 2 (P 2 ) /^tanh(M) 



(to - V) 2 

?; e(M-|A 2 (p 2 )|) | e(HHAi(p^)l)e(|A 2 (p^)|-|^ + i1) '| (36) 

/^_|A 2 (p 2 )|2 .^-|Ai(pi)l 2 x /\A 2 (p 2 )\ 2 -(co + Vy ] 

u + V 
IT 



oo 

I 2 {V, pi) = A 1 (p 1 )A 2 (p 2 ) J duo (tanh (^\ - tanh 



K sgn(uj)sgn(uj + V)B - |Ai(pi)]) 9 (\u + V\- |A 2 (p 2 )|) 

/c;2-|Ai(pi)| 2 J(uj + V) 2 -\A 2 (p 2 )\ 2 ' 



WM= I du ( tanh - tanh (JL)) Ml" ~ V \ @ M Z ^ffiOD x 



9 (|a;-\/|-|A 2 (p 2 )|) 
W -\/) 2 -|A 2 (p 2 )P 



(38) 



WPi) 



^tanh ( JL) {uJ - V )\.\ (ptttax 



w e(M-|A 2 (p 2 )|) e (H - |Ai(pi)|) e (|A 2 ( P2 )| - \u - v\) 

X = — h 



uj 2 - |A 2 (p 2 ) 



|A a (pO| 2 J|A 2 (p 2 )| 2 -(c-^ 



(39) 



The singularities of functions I m (V, pi) appear after the integration over the frequency 
only if two square roots (multiplied by each other in the denominators of the integrands) are 
equal to zero simultaneously. It is possible for certain values of frequency and voltage and the 
singularities turn out to be located at momentum-dependent points \V\ = ||A 2 | ± |Ai||. It 
follows from Eqs.(|36|)-(|3"9"|) that the expressions for singular parts of I m (V, pi) are as follows 



I x = -^/\A l A 2 \sgn(A 1 A 2 ) 



tanh 



|Ai| 
2T 



+ tanh 



' [A 2 [ 
2T 



x 



x In 



■|Ai| -|A 2 ||-tt 



tanh I 



'JAil 
2T 



tanh | 



|A 2 | 
2T 



e(|V|-||A 2 |-|Ai 



(40) 
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h = -^\A 1 A 2 \sgn{A 1 A 2 )sgn{V) i sgn (|Ai| - |A 2 | 



— tanh 



tanh 



(All 
2T 



A l 


+ 


V 


sgn (|A 2 






lnHVI - ||A 2 | - |Ai||| 


2T I 



+7T 



tanh 



AJN , /|A 2 
+ tanh 



2T 



2T 



e(W\ 



h = - V|AiA 2 |s^n(V) (| Ai| - |A 2 |) 



tanh 



[Aij + |V^n(|A 2 |-|Ail) ' 
2T 



In V 




(41) 



— 7T 
1 

~2 

+ 7T 



tanh 



Ai|\ , /IA21 
+ tanh 



2T 



AiA 2 



tanh 



tanh 



\A X \ 
2T 



tanh 




e(|y|-|A 1 |-|A 2 |) 



(42) 



+ tanh 



|A 2 | 
2T 



In 



IAJ - IA2II + 



e(|F|-||A 2 |-|Ai| 



(43) 



According to Eqs.(|6|), (0), for obtaining the current- voltage relation in the vicinities 
of singular points one should carry out the integration of expressions (HDp-(I3) over the 
momentum directions and consider the correspondent expressions for j m (V) . It follows from 
this integration that nonanalytical behavior of the I-V curves takes place only for the values 
\V\ in the close vicinities of the extremal points of the expressions ||A 2 (p 2 )| ± |A 1 (p 1 )||, 
which are considered, for example, as functions of pi. The characteristic behaviors of the 
I-V curves near these extrema turn out to be strongly dependent upon the type of the 
extremal point. Comparing the expressions one fi n ds the same singular behaviors 

for pairs of functions ji(V), ji(V) and jiiy) , JziV) (disregarding the differences in signs 
for the moment). Due to this fact it is sufficient to describe, for example, only the singular 
points for ji(V) and jz(V). 

It is convenient further to examine the singular points for the conductance G = dj x /dV. 
It can be shown, that only the derivatives of 0- and logarithmic functions with respect to the 
voltage are associated with nonanalytical behavior of G(V). In the former case one obtains 
5-function, which in fact reduces the integration over the Fermi-surface to the integration 
over the line on the surface. Corresponding terms Gi^iV) of functions Gi^iV) = dji^/dV 
may be represented as follows: 



G 1 = -sgn{V) dl 



K- 



Pi 



|Vp X (|Ai| - |A 2 | 



G, 



dl 



K + ( Pl ) ^n(AiA 2 ) 



pi 



lAJ + IAa 



(44) 
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Here / is the local coordinate along the line \V\ = ||Ai(pi)| ± | (P2) 1 1 on the Fermi 
surface (sign plus corresponds to G 3 and minus - to G\). The functions K ± are determined 
as follows: 



tanh ± tanh (§f) | yj\ A 1 A 2 \sgn{A 1 A 2 )^D 



K* (pO = 7T v ^ 7 % : : ^— (45) 

lPU 2(2tt) 3 1 ' 

We consider below various types of extrema and get respective behaviors of G m (V). Let 
the function 1 1 A2 1 — |Ai|| take maximal or minimal value at the point pi = po on the Fermi 
surface and in the vicinity of this point has the form 

||A 2 | - |Ai|| = a ± (bp\ + cpl),a,b,c> . (46) 

Here pi, p 2 are the local orthogonal coordinates in the vicinity of the point po- Since the 
function ||A 2 | — |Ai|| comes in the expressions for j 14 as the argument of O-function and 
in the formulae for j 23 as an argument of the logarithmic function, we obtain two different 
singular behaviors near the value |V| = a: 

5G l \\ V \=a= s^(AxA 2 ) | Pl=P0 5G 4 \\ V \= a = T^=K~ (p ) (47) 

ybc 



G 2 = sgn(A 1 A 2 ) | Pl=Po G 3 = ±4=^" (Po) ^ \\V\ - a\ (48) 

ybc 



The notation for a jump of the conductance 8G ||y|= a = <j(|V| > a)— G(\V\ < a) is introduced 
here. 

The singularities coming from the maximal and minimal values of the quantity | A 2 1 + 1 Ai | , 
when one has |A 2 | + |Ai| = a ± (bp\ + cp 2 ), a,b,c> 0, are described analogously: 

$G 2 \\ V \=a= -s^n(AiA 2 ) | P1=P0 5G 3 \\ V \=a= ^^ K+ (Po) ( 49 ) 

Gi = -s<7n(A!A 2 ) |. i= . o G 4 = ±^K + (po) sgn(V) In \\V\ - a\ (50) 

In the case of a saddle point of the function ||A 2 | — |Ai|| one has near this point 

||A 2 | - |Ai|| = a + bpl - cpj, a, b,c>0. (51) 
The respective singular parts of the conductance read 

G 1 = s^n(A 1 A 2 ) | P1=P0 G 4 = -±=K~ (p ) sgn(V) In | \V\ - a\ (52) 

ybc 



G 2 = sgn(A 1 A 2 ) | Pl=Po G 3 = -^=sgn(\V\ - a)K~ (p ) In 2 \\V\ - a\ (53) 

ny be 
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Analogously, for the saddle point 



| A 2 | + | Ai| = a + bp\ — cp 2 , a, b, c > 



(54) 



we get 



G 2 = -sgniA^) 



IPl=P0 



G 3 



K + (p ) \n\a-\V\ 



(55) 



Gl = -sgn(A 1 A 2 ) | Pl = po G A = ^n(V)s g n(\V\ - a) R+ ^ ||y| _ q| 

7TV0C 

The jumps in conductance described by Eqs.Q4"7|),(f^) correspond, of course, to the kinks 
on the I-V curves. Logarithmic divergences in conductance described by Eqs.(fi8|)X|50"|),(|52|) 
and (KBp, result in step-like points on the I-V curves (note that on both sides of these 
singular points G has the same sign). At last the terms containing logarithm squared in 
Eqs.([53|), (|56f) describe the cusps (beak-like points) on the current- voltage characteristics. 



They appear in the case of saddle points of the functions ||A 2 | ± |Ai|| after the integration 
of logarithmic singularities for I m . 

Let now the quantities | |A 2 | ± | Ai| | have the extremal values on some line I on the Fermi- 
surface, rather than on isolated points as it was suggested above. Then, for example, in the 
vicinities of maximal or minimal values of ||A 2 | — |Ai|| one has 

||A 2 | - |Ai|| = a±bp\ , a,6>0, (57) 

where pi is the local coordinate on the Fermi- surface orthogonal to the extremal line 

In this case the inverse square root singularities in the conductance appear from one side 
of the voltage value \V\ = a: 

Gl = sgn(A 1 A 2 ) \ t G 4 = - ^2 6 (± (\V\-a)) f dl ^= , (58) 



C , = , t »(A,A ! )|,^. T e( f | -° )1 / dl^. (59) 

\ \\V\ - a\ IIA . i A . VO 

Vxl>0 

These singularities correspond to the vertical slope of the I-V curve from one side of the 
voltage value |V| = a for each of four terms presented in the total expression for the tunnel 
current. For example, in the case of maximum on the line I there are the vertical slope of 
the curves ji(V) and j&(V) at \V\ = a from side |V| < a, and of the curves j 2 (V) , jziV) 
from side \V\ > a. 

Analogously, in the case of maximal or minimal value of the quantity |Ai| + |A 2 | on the 
line I, when near this line one has 

|Ai| + |A 2 | = a±bp\ , a,6>0, (60) 
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the singular behavior of the conductance is described as follows: 

G2 = - sgn( A 1 A 2 )| rG 3= e( t (|1/| " a)) / (61) 

' I VI - a] J 2Vo 



G? 1 = -^n(A 1 A 2 )| r G 4 = T-^J=e(T(|V|-fl)) / (62) 

^/MVI — «| |Ail+ ^ 2l=a Vb 

Vxl>0 

We have considered above different superconductors or at least different crystal orienta- 
tions from both sides of the junction, when not only sum |Aj| + |A 2 | may have extremal 
points or lines on the Fermi surface, but the difference |Ai| — |A 2 | also depend upon the 
momentum direction and may have extrema. As a result one can't take directly the limit 
Ai = ±A 2 = A in the expressions written above for singularities, which are associated 
with extrema of the difference |Ai| — |A 2 |. For this particular case quantities have not 
singular point at V = |Ai| — |A 2 | = 0, while for G 2; 3 in the vicinity V = ( more exactly, 
for \V\ < T) one gets instead of Eqs. (H),(0),(||)' 

G -±C -± ln|V| f d2SlVxl D(u) |A| f6tf 

° 2 ~ ±Gs ~ ± ^ J J^f^ (Pl) • (63) 

For isotropic s-wave superconductors this term is exponentially small at law temperatures 
T < A. In contrast to this case, for anisotropically paired superconductors the expression 
in Eq . (|63|) manifests power law temperature behavior. For instance, for the line of nodes of 
the order parameter ( when in its vicinity |A(p)| = 6|pi| ) it follows from Eq.(|63"D at low 
temperatures G 2j 3 oc Tin |V|. 

The results of numerical calculations for j m {v) (v = V/A Q ) for the case when there is 
no surface pair breaking from both sides of the tunnel barrier are shown in Figs. 2, 3. In 
Figure 2 the junction between isotropically and anisotropically paired superconductors is 
considered: Ai = A o cos(20) , A 2 = A /2 = const. Here - is the azimuth angle in xy - 
plane of a tetragonal superconductor (axis z is parallel to the boundary plane). The Fermi 
surface is assumed to be cylindrical for anisotropically paired superconductor (with d-wave 
pairing). We let for the barrier transparency D oc cos 2 (0) and Aq/(2T) = 0.5. In this 
case the singular points on the I-V curves are only the points of maximum values of the 
quantities ||Aj(0)| ± |A 2 ||, as the minimum corresponds to the zero value of Ai. Note that 
for A 2 = 2Ao and for the same Ai(0) only minima of the quantities 1 1 Ai(0) | — | A 2 | | would be 
of importance. In Figure 3 the case of two identical anisotropically paired superconductors 
is described: Ai = A 2 = A o cos(20). All functions j m (f) are normalized to the value |ji(0)|. 

B. Crystal Orientations with Surface Pair Breaking from One Side of the Barrier 

Plane 

Let us discuss now the tunnel junction between two anisotropically paired superconduc- 
tors, considering a gradual change of the crystalline orientation of one of them relative to 
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the barrier plane and retaining the condition A 2 (p) = A 2 (p) to be fulfilled only for the 
second superconductor. According to Eqs. fl34|) , (|35|) , there are no square root divergences of 
the propagators for the first superconductor, taken on the barrier plane for an intermediate 
crystal orientation. Hence, the singularities, which were found in the previous subsection, 
must begin to become smooth and subsequently, for large enough deviations from the initial 
orientation, will disappear. At the same time, as it was already mentioned above, some new 
characteristic singular points on the I-V curves appear in this case. One kind of them turns 
out to be associated with the existence of regions on the Fermi surface with the opposite 
signs of the order parameter A l oo (p). Then the zero energy quasiparticle bound state oc- 
curs at the boundary plane. Other singularities appear if the additional quasiparticle bound 
states localized near the barrier plane appear due to the particular form of the spatial depen- 
dence of the order parameter ||. As it is seen from Eqs.([34]), (|35|), the terms in propagators 
containing the factor 1/uo are of importance for the former case. For the latter case the 
contribution to the current from additional poles of the propagators has to be taken into 
account. So, we let the singular parts of the propagators for the first superconductor be 
represented by Eqs.fl34"|), (pop, and for the second superconductor we use Eq. (|33|) . 

It is very essential that in the case when zero energy bound state takes place only from 
one side of the junction (while from the other side there is no surface pair breaking), the 
respective singular contribution comes only to the quasiparticle current, not to the Josephson 
effect. It can be shown under the same conditions as are supposed to be fulfilled above. The 
point is that for the singular parts of the propagators, associated with the mid-gap state, 
one has the following property f s (f>) = / s + (p)- Opposite to this equality, the relation 
f(p) = — / + (p) holds for the orientations when there is no surface pair breaking. Due to 
this reason the corresponding singular parts of ji, j 2 reduce to zero. It is not the case for 
the bound states with nonzero energy, since the function Qf{p) is complex in contrast to 
the real function B f (p) (Bj = B) = B f , Q + f = Q* f ). 

Substituting the expressions (0), (P5|), ( |55D into Eqs.(|TUp, ([□]), we get for / 3i 4 the 
singular points which appear after the integration over u: 



( v 

h(V, pi) = -B gl {pi) tanh — 



Qgi(Pi)y\^2\sgn(V) 
V2 



V2T 



|y|e(|v|-|A 2 (p 2 ) 



(tanh tanh J 



|A 2 (p 2 )P 

e (\v\ - h - \a 2 \) 



|V|-/ii-|A 2 | 



+ 



tanh 



|A 2 | 
2T 



tanh 



'V 

2T 



|V|-||A 2 |-/ii|)s</n(|A 2 |-/H) 



(64) 



l4(V,pi) = B fl i(pi)tanh 



\v\\ |v|e(|A 2 (p 2 )| - 



Q g i(pi)y|A 

V2 



2T 



l (tanh ( 



2T \ 2T 
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The nonanalytical terms of square root type presented in Eqs.(^), fl55|), are omitted in 
Eqs.(|64]), (|65|), since in this case they result in jumps or divergences only for the derivatives 
of the conductance, not for the current or for the conductance itself. 

Further integration over the Fermi surface may lead to different kinds of the singular 
points on the I-V curve, since there are various possibilities for the behaviors of the order 
parameter A 2 and function hi with the momentum and types of the corresponding extrema 
on the Fermi surface. We consider firstly the singular points of the I-V curves coming from 
the terms of l/tu-form in the expressions for the propagators of the first superconductor 
and confine ourselves by two examples. In the particular case of isotropic s-wave second 
superconductor there are the inverse square root singular point on the current- voltage char- 
acteristic: 

]3 = - tanh(-j _ J DB tl ( Pl ) , (66) 

V ' ' v x i>0 



. f\v\\ \v\e (\a 2 \ - \v\) f d 2 s lVxl 

j 4 = tanh — / — — — DB gl Pl . 67 

So, the quantity j'3 diverges at |V| = A 2 from the side \V\ > A 2 , and the quantity j 4 - 
from the side \V\ < A 2 . In the case A 2 = (S-N junction) we get from Eq.(|66|) that the 

conductance G3 oc (T cosh 2 (V/2T)^j . Then for small voltage \V\ <C T it follows G3 oc 1/T 
and we obtain the zero-bias anomaly of the conductance at low enough temperatures (see also 



21|). The divergence of G3 in the zero-temperature limit takes place only in the idealized 
system without taking account of the factors, which broaden delta-peaks in the quasiparticle 
density of states. 

Further, if the anisotropic order parameter A 2 has an extremal line on the Fermi surface 
and in its vicinity one gets 

|A 2 | = a±bpj, a,b > , (68) 

then the kind of the singular point depends on the behavior of B g i(pi) near this line. If this 
function doesn't vanish in the vicinity of the extremal line, the logarithmic divergences take 
place for j 3 for the line of maxima and for j 4 for the line of minima: 

J 3 ,J4 oc y^tanh (j^j In ||V| - a\ . (69) 

In the case of vanishing of the function B g i(pi) on the line, when in its vicinity this 
function has the form 

B g i(pi) = P\pi\ , (70) 



18 



one-sided vertical slope for functions j^i appears at \V\ = a: 

( a \ 9 (\V\ - a) „ \fa . , / a \ (-) ((/ 



Gs oc tanh . 

3 6 V2T 



, G 4 oc — — tanh , 

b \2T 



\V\) 



\v\ 



(71) 



Now let us examine the singular points coming from the terms describing the additional 
pole in the propagators for the first superconductor. In contrast to the contribution from 
the midgap states, even in the case of surface pair breaking from one side of the barrier 
plane the quasiparticle states bound to the tunnel barier with nonzero energy contribute 
not only to the quasiparticle current but also to the Josephson effect. It takes place due to 
the nonzero values of ImQf(p), while ImBf(p) = 0. Positions of respective singular points 
on the I-V curves turn out to be associated with the extremal points on the Fermi surface 
of quantities hi ± IA2I. Considering, for example, the line of extrema I for the quantity 



hi + |A 2 | =a±bp\ 



a,b> 



we find the logarithmic divergences and jumps for j T , 



(72) 



(73) 



Here the upper (down) line in all four expressions corresponds to the upper (down) sign 
in Eq.([F2|). We introduce the notations: 



7 7 -M + i ln||y| " a| 7 7 - M + san(V) I ^ " 

3i, ^- M f, 9 \nQ{\V\-a) ' 32,3s- M f ^sgn{V)< ln ||y|_ a | 



Ml 



1 



16^7^ 



^i>0 



I A, 



-Qgl(Pl) 



tanh (*lU tanh 

\2T \ 2T 



(74) 



Mf= fdl^D JM 8 ^) ImQ fa 
f J v n V b 16V27T 3 f 



If the function ||A 2 | — hi\ has the line of extrema 

IIA2I — hi\ — a ± bp\, a,b>0 , 

then one gets 

ji A = Mj >g { T Q (± (|A 2 | - hi)) \n\\V\ -a\ + tt0( T (|A 2 | - hi))Q(a - \V\)} 



(75) 



(76) 



(77) 



j 2 ,3 = M^n(y){e( T (|A 2 |-/i 1 ))lii||y|-a|±7re(±(|A 2 |-/i 1 ))e(a-|y|)}. 



(78) 



So, for positive (negative) value of |A 2 | — hi the function j 2 3 (ji^) has logarithmic 
divergence only for the line of maxima of ||A 2 | — hi\, while j'14 ^2,3) has such divergence 
only for the line of minima of the same quantity. 
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As it was mentioned above, the function h(f>) may manifest the nonanalytical behavior, 
for example, at the momentum direction, for which the bound state at the boundary disap- 
pears. Taking this possibility into account, we suppose now that the function hi + |A 2 | has 
near some line on the Fermi surface the following nonanlytical form 



hi + |A 2 | = a + (66(p) + cQ(-p))p ,a > . 



(79) 



Particular cases b (or c) — > oo just correspond to the absence of the bound state for p > 
(or p < 0) on this line. Then one-sided vertical slope for functions j m appears at \V\ = a: 



_ Q(a-\V\) 

Ui, — (jt4 — 

y/"-\V\ 

We introduce the notations: 

V X 1 



sgn(V)P+ g , G 2 , -G 3 



e( iv|-q) 

|V|-a 



P 



f,g 



r// ( I - - 



v x i>0 



cj v n 



DJ\A, 



| Qgl(Pl> 



, , /|A 2 | 

tanh — — - ± tanh — — 
\2T/ I 2T 



(80) 



(81) 



P ± 



Wrl>0 



CJ Vfi 



-sgn(A 



2| ^7I^ /mg/l(lJl 



tanh [ — i- 1 ±tanh f J—^- 
l2T/ I 2T 



(82) 



If the function \hi — |A 2 || has the form (|79|) near some line, then also one-sided vertical 
slope for functions j m appears at |V| = a: 



G M = 55fn(V A )s5fn(|A 2 |-/ii) 



e((a-|\/|)(|A 2 |-^)) 
a - |y|)s^n(|A 2 | -/n) 



P 



(83) 



G 2 ,3 - ^(ri-a)^-^) 

The Eqs.(|64|), (|65|) are not appropriate for direct consideration of the case A 2 = for 
a finite voltage, since the voltage is supposed to lie near the corresponding singularity for 
each pole-like term. For S-N junction it follows from Eqs. (0), (|10|), fllll) , that in the case 
of extremal point of the quantity hi {hi = a ± bpl, a, b > ) there is the low-temperature 
anomaly in the conductance for \V\ = a: G 3)4 oc 1/VUT. 

The results of numerical calculations for j m {v) (v = V/A ) for the case when there is 
surface pair breaking from one side of the tunnel barrier are shown in Figs. 4.1, 4.2. The 
tunnel junction between <i-wave and isotropic s-wave superconductors is considered under the 
conditions: Aloo = A cos(20 — 2(f> ), A 2 = 0.2A = const. For the rf-wave superconductor 
we choose the same parameters as earlier (see Fig. 1): (f> = n/9, T = 0.45T c i, A /(2T) = 2. 
Here T cl is the critical temperature for the superconductor with d-wave pairing. We let for 
the barrier transparency D oc cos 2 (</>). 

There are the inverse square root singular points on the curves j'3, j'4 for V = A 2 = 0.2A . 
For V = h m — A 2 = 0.5A there are logarithmic divergences of ji, j 4 and jumps of j 2 , j 3 . 
For V = h m + A 2 = 0.9A there are logarithmic divergences of j 2 , j 3 and jumps of ji, j 4 . 
At voltages V = h ed ± A 2 = 0.83A , 0.43A there are kinks on the I-V curves, though some 
of them are feebly marked. All functions j m (v) are normalized to the value ji(0). 
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C. Crystal Orientations with Surface Pair Breaking for both Superconductors 



Let now the both superconductors have the intermediate orientations relative to the 
boundary. Then one should make use the expressions fl34|), (|35|) for the singular parts of the 
propagators on both banks of the barrier plane. Substituting these expressions into Eqs. ([§])- 
(|ll|), we obtain several kinds of the singularities. They come from the multiplications of the 
terms with poles by each other and from the products of the pole-like term and the square 
root nonanalytical term, represented in Eqs.(^), fl35|). All these terms result, for example, 
in the following singular part of ^3 

h = ^ tanh Q^i (pO^Oi (Pa) v^-Aic (p 2 )e^| - |A 2oo (p 2 ) |) + 



imc 2 (ft,) ^/al, ( P2 ) - v*e (| a 2oo (p 2 ) | - - |g 92 (p 2 ) s (\v\ - h 2 (p 2 ))+ 

+ (P 2 - P 2 )) + l{ (tanh (^M±Z) - tanh {^M\\ Q gl ( Pl ) . 



{-\Q 9 2 (P2) 5 (|Ai (Pi)+ V\-h 2 (p 2 )) + i?eC 2 (pa) v 7 ^ (pi)+^) 2 -AL (p 2 



6 (\h ( Pl ) + F| - |A 2oo (p 2 )|) + /mC 2 ( P2 ) v Aloe (P2) - (h (pi) + V) • 

e (|A 2oo (p 2 )| - |/n ( Pl ) + y|) - (v ^ -v)) + (p 2 ^ p 2 )} + (1^2) , (85) 

Analogues singularities appear in the expressions for / 12 ,4- 

Here we have omitted the joint contribution from the zero frequency singularities of the 
propagators from both sides of the junction. For finite voltages this contribution is equal to 
zero in the idealized system in question, though it may be of importance for low voltages 
less or comparable with the width of the broaden delta-peak in the density of states. For dc 
Josephson effect this term leads to the low temperature anomaly in the critical current and 



to the possibility of — n phase transition I 22|j23 |. 



Further integration over the momentum direction is performing below for the particular 
cases of the lines of extrema for the order parameter | A 2oo (p 2 ) | or for quantities \h% ± | A 2oo | |, 
h,2, |/ii±/i 2 |. Let , for example, the order parameter |A 2oo (p 2 ) | or the quantity \hi ± |A 2oo || 
have the line of extrema of the form 

|A 2oo | , \hi ± |A 2oo || = a ± bp 2 , a,b>0 , (86) 

and the quantity £>i(pi), or |A 2oo |, Qi(pi) is not equal to zero on this line. Then the terms 
in the conductance £a j2j 4 acquire the following logarithmic singularities 

sgn(V)Gi,G 2 ,sgn(V)G 4c oc y^rn||V| - a\ . (87) 

At the same time the logarithmic divergence (0) of G3 appears only in the case of the line 
of maxima, due to the relation ImC = under the condition \u\ < min (|Aoo(p)|, |Aoo(p)|). 
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For the particular orientation, when one has A 2oo (p 2 ) = — A 2oo (p2), the logarithmic diver- 
gence of G4 remains only for the line of minima. Note, that the square root behavior of 



the propagators of the form y A^ (p) — uj 2 result here in the logarithmic divergences of the 
conductance. For some other cases these terms may result in the kinks for a conductance 
or in the divergences of its derivatives. We do not consider these kinds of singularities here. 
If there is the line of extrema of h 2 

h 2 = a±bp 2 , a, b > , (88) 

then the following inverse square root divergences appear on the I-V curves: 

h, TH = T^tanh f^) 9 ~n == =^ / dl^B^^ReQ ^2)^= , (89) 
16tt 3 \2TJ J\\ V \-a\ J >0 v f i Vb 



32, -h = ^tanh (¥-) 9( ^ | ( 'y l ° )} f dl^B^i^ReQ^)^ . (90) 
V2T7 ^\\V\-a\ J >0 v fl Vb 

In this section coefficients Bf 2 , Q*f 2 originate from the corresponding expression for / + 
for second superconductor, while Bfi, Qf 1 come from the expression for / for the first 
superconductor. 

In the case of nonanalytical behavior of the function h 2 near some line one has, for 
example, 

h 2 = a + (66(p) + cQ{-p))p , a > . (91) 
Then the logarithmic divergences of occur at |V| = a: 



3\, M = ln|l ^4 a ' tanh (^) / dl^DBf^ReQf^) Q - \) , (92) 

^ 'v x i>0 fl 

and the following jumps appear in j 2 ,3 at |V| = a: 

-j 3 =^ tanh Q J ^fDB^MReQ^M ( ° & = °» 6 » + 

^!>o /1 VII 

+ e » a -' y " C ») . (93) 

At last, let us consider the case of the line of extrema of the quantity \h\ — h 2 \: 

I hi - h 2 \ = a ± bp 2 , a, b > . (94) 
Then one gets divergences on the I-V curves of the inverse square root types: 

T&(±(a-\V\)) r dl Ar 
167ry||y|-a| J >0 Vb 
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■ ■ sgn(V)Q(±(\V\-a)) f dl 

•72> = 7= / -7=Nf, g , (96) 



87T 3 »/ijvi - oj „ J Vb 

V 1 1 1 1 «S!l>0 

where we introduce 

Nf, 9 = V f-DRe (Q f , gl Q* lg2 




(97) 



In the case when the quantity \hi — h 2 \ near some line behaves like ([JI]), the logarithmic 
divergences of ji j4 appear at \V\ = a: 

v xl >0 

and the following jumps take place in j 2 ,3 at \V\ = a: 

■ ■ S ^ n{ y) [ mm ( Q((\V\-a)b) , Q((a-\V\)c) \ 

j2 '- j3= ^^ y — \b\ — + — m — ) (99) 

The results for the line of extrema of the quantity hi + /i2 follow from fl9~4])-(|9"PD after the 
substitution h 2 — > —h 2 , ji — > —ji, Q*f 2 — > Q/2- One should note, that in addition to the 
presented in this subsection singular points there are the analogous ones which may be 
described by the same equations (|86|)-(p3|) after the interchange 1 <-> 2. 

The results of numerical calculations for j m (v) (v = V/Aq) for the case when there is 
surface pair breaking from both sides of the tunnel barrier are shown in Figs. 5.1, 5.2. 
The tunnel junction between two identical c?-wave superconductors is considered for the 
particular case of "mirror" junction, when the barrier is a reflection-symmetry plane of the 
superconducting electrodes: A loo (p!) = A 2oo (pi) = A cos(20 — 20 o ). As earlier we let (see 
Fig. 1) (fi = 7r/9, T = 0.45T C , A /(2T) = 2, D oc cos 2 (0). There are the inverse square root 
singular points on the curves ji, 2,3,4 at voltage V = 2h m = 1.4A . For V = 2h e d = 1.26A 
there are logarithmic divergences of ji, j'4 and jumps of j 2 , j'3. At voltage V = 0.635Ao 
curves ji, j 2 have kink-like behavior, which may be associated with the interplay between 
the contributions of mid-gap states from one side of the junction and of continuous quasi- 
particle spectrum from another side: the function min (| cos(20 — 2</> )|, | cos(20 + 20 o )|) has 
maximal value 0.635 for the direction = ±0.79 at which Bf ^ 0. 



V. CONCLUSIONS 

As we have shown above, there may be a large variety of different types of nonanalyti- 
cal points on the I-V curves for tunnel junctions in anisotropically paired superconductors. 
The singular behavior differ essentially from that which is characteristic for case of s-wave 
isotropic superconductors. The quality of the barrier plane may have a great influence on 
the manifestations of the singular points in question. For example, instead of the divergences 
of the current, which are obtained for the idealized system, there are peaks, whose practical 
finite magnitudes are sensitive to the values of elastic and inelastic scattering processes, 
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roughness of the barrier plane, and even the finite value of the junction transparency. The 
inelastic scattering processes are common for removing the singularities of the tunnel current 
in s-wave isotropic and in anisotropically paired superconductors. The elastic scattering pro- 
cesses are pair breaking just for the anisotropically paired superconductors. Besides, these 
factors together with the quality of the barrier plane and the finite value of transparency 
may broaden the delta-peaks in quasiparticle density of states and, subsequently, wash out 
the corresponding peaks in the tunnel current. Nevertheless, the characteristic behavior 
of the I-V curves considered above is observable under the certain realistic conditions and 
may be employed as a sensitive test for identifying the anisotropic types of pairing in the 
superconductors, in particular, different signs of the order parameter on the Fermi surface. 



ACKNOWLEDGEMENTS 

One of us (Yu.S.B.) is grateful to D. Rainer and J. A. Sauls for stimulating discussions and 
for sending the preprints of their works before publication. This work was supported by the 
grant No. 96-02-16249 of the Russian Foundation for Basic Research. A.A.S. acknowledges 
Forschungszentrum Jiilich for financial support. 



24 



REFERENCES 



[1] A.I. Larkin, and Yu.N. Ovchinnikov : Zh. Eksp. Teor. Fiz. 51, 1535, (1966). 
[2] N.R. Werthamer, Phys. Rev. 147, 255 (1966). 

[3] A. Barone, and G. Paterno, Physics and Applications of the Josephson Effect (Wiley & 
Sons, N.Y.) (1982). 

[4] V.B. Geshkenbein, LA. Larkin, and A. Barone, Phys. Rev. B 36, 235 (1987). 
[5] M. Sigrist, and T.M. Rice, J. Phys. Soc. Jpn. 61, 4283 (1992). 
[6] A. Millis, D. Rainer, J.A. Sauls, Phys. Rev. B 38, 4504 (1988). 
[7] Yu.S. Barash, A.V. Galaktionov, and A.D. Zaikin, Phys. Rev. B 52, 665 (1995). 
[8] L.J. Buchholtz, M. Palumbo, D. Rainer, and J.A. Sauls, J. Low Temp. Phys. 101, 1099 
(1995). 

[9] M. Sigrist, and T.M. Rice, Rev. Mod. Phys. 67, 503, (1995). 
[10] D.J. Van Harlingen, Rev. Mod. Phys. 67, 515, (1995). 
[11] Ch. Renner, and 0. Fischer, Phys. Rev. B 51, 9208, (1995). 

[12] Yu.S. Barash, A.V. Galaktionov, and A.D. Zaikin, Phys. Rev. Let. 75, 1675 (1995). 

[13] E. Riedel, Z. Naturforsch. 19A, 1634 (1964). 

[14] C.-R. Hu, Phys. Rev. Letters 72, 1526 (1994). 

[15] J. Yang, C.-R. Hu, Phys. Rev. B 50, 16766 (1994). 

[16] M. Matsumoto, H. Shiba, Journ. Phys. Soc. Jpn. 64, 4867 (1995). 

[17] Y. Nagato, and K. Nagai, Phys. Rev. B 51, 16254 (1995). 

[18] A.V. Zaitsev, Zh. Eksp. Teor. Fiz. 86, 1742 (1984) [Sov. Phys. JETP 59, 1015]. 
[19] L.J. Buchholtz, M. Palumbo, D. Rainer, and J.A. Sauls, J. Low Temp. Phys. 101, 1079 
(1995). 

[20] M. Matsumoto, H. Shiba, Journ. Phys. Soc. Jpn. 64, 3384 (1995). 
[21] Y. Tanaka, and S. Kashiwaya, Phys. Rev. Letters 74, 3451 (1995). 
[22] Yu.S. Barash, H. Burkhardt, D. Rainer to appear in Phys. Rev. Letters. 
[23] Y. Tanaka, and S. Kashiwaya, Phys. Rev. B 53, R11957 (1996). 



25 



o,a 




angle <|> (radian) 
Figure 1 




0,0 0,5 1,0 1,5 2,0 2,5 

voltage v 



Figure 2 




.3 i — ' — i — ' — i — ' — i — ' — i — ' — i — ' — i 

0,0 0,5 1,0 1,5 2,0 2,5 3,0 

voltage v 
Figure 3 




0,0 0,2 0,4 0,6 0,8 1,0 1,2 

voltage v 
Figure 4.1 



CD 



O 




i-0,5- 



-1,0- 



-1,5 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 1 

0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 




O0,0 



-0,5 - 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0,0 0,2 0,4 0,6 0,8 1,0 1,2 



i 1 1 1 1 1 
1,4 1,6 1,8 



voltage v 
Figure 5.1 



10 



8- 




-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 

voltage v 
Figure 5.2 



